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Abstract 

We update the theoretical analysis of the CP-violating asymmetry in the decay Kl — ► 
7r + 7r~e + e~, relying on chiral perturbation theory and on the most recent phenomenologi- 
cal information. With the experimentally determined magnetic amplitude and branching 
ratio as input, the asymmetry can be calculated with good accuracy. The theoretical 
interpretation of the sign of the asymmetry is discussed. 
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1. A large CP- violating asymmetry in the decay Kl — > 7r + 7r~e + e~ was originally pre- 
dicted by Sehgal and Wanninger |T|. The effect is almost entirely due to indirect CP 
violation in K° — K° mixing. The predicted asymmetry in the angle between the 7r + 7r~ 
and the e + e~ planes has been confirmed experimentally || [J| . The purpose of this note 
is to update the theoretical analysis of the asymmetry |1|, f|, ^ ||, [7], || using the most re- 
cent phenomenological input available and employing the methods of chiral perturbation 
theory. The latter will mainly be invoked to estimate higher-order corrections but also to 
interpret the observed sign of the asymmetry. 

2. The amplitude for the decay Kl{j>) — > it + (pi)it~ (p 2 )e + {k + )e~ (kJ) is expressed in 
terms of three invariant form factors E\,E2 (electric) and M (magnetic): 

A(K L -> Tr+TTe+O = ^n(k-)'fv(k+) V, 

q 2 

V„ = iE lPlll + iE 2 p 2ll + Me fll/pa p u 1 p p 2 q a (1) 
q = k + + fc_ . 

We use the set of variables originally introduced by Cabibbo and Maksymovicz || for 
K e 4 decays. With this choice, the form factors Ei,E 2 ,M depend on s n (invariant mass 
squared of the two pions), q 2 (invariant mass squared of the leptons) and 9 n (angle of 
the 7r + in the n + 7i~ center-of-mass system with respect to the dipion line of flight in the 
kaon rest frame). The two remaining Dalitz variables are 8i, the corresponding angle for 
the positron, and $, the angle between the dipion and dilepton planes in the kaon rest 
frame. 

After integration over four of these variables, the differential decay rate with respect 
to $ assumes the general formQ ]10], |TT 



— = h + I 2 cos $ + I 3 sin $ + J 4 cos 2$ + J 5 sin 2$ . (2) 
Under a CP transformation 

cos $ — ► cos $ 
sin $ — > - sin $ (3) 

so that non-zero J3, J5 signify CP violation. It turns out that ^3 is very small in the 
standard model, being sensitive to direct CP violation only [|J. The quantity of interest 
here is 1$ that is almost exclusively due to indirect CP violation |3j]. A convenient measure 
of this term is the asymmetry 

4/5 

A C p = (sgn(sin $ cos $ = — _ . (4) 

T(K L — > 7r + 7r e+e ) 



L We use the conventions of Ref. [fl3j| , in particular £0123 = 1 and the metric (H ) 
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In terms of the invariant form factors denned in (JJ), the asymmetry is given by (T^, [TT 



n'2 



AcP 1287T 8 M K oT(K L -> 7r+7r-e+e-) ^ 



/ dupse^rfpZKk* {Pl P2) ' [ K ^ sgn(sin $ cos $)Im[(E 1 - E 2 )M* 
J q 



where dups is the invariant phase space integration measure. We will also consider the 
asymmetry for certain cuts (in q 2 ). In this case, both the phase space integration in 
the numerator of @ and T(Kl — > 7r + 7r~e + e~) in the denominator must be modified 
accordingly. 

When CP is conserved we have E 2 (pi,p 2 ,q) = Ei(p 2 ,Pi, q)- Therefore, only CP- 
violating contributions to E\, E 2 matter for the numerator in (|5|). 

The most recent published result for the asymmetry comes from the KTeV- 
Collaboration 0: 

A C p = 13.6 ± 2.5 (stat) ± 1.2 (syst) % . (6) 

The preliminary result from NA48 || is fully compatible with this value. 

Before delving into the theoretical analysis, we comment on the observed sign of the 
asymmetry @. With our conventions |TTJ and with the assumption that lm[(Ei — E 2 ) M*] 



has a unique sign all over phase space, the theoretical expression (Jg) implies 

sgn A CP = sgn lm[{Ex - E 2 )M*] . (7) 

As we shall show in the sequel, Im[(£/i — E 2 )M*} does in fact have a unique sign and 
must therefore be positive according to Eqs. (|G[) and (0). 

3. At lowest order in the low-energy expansion, the electric amplitudes correspond to 
pure Bremsstrahlung: 



-2ieA tree (K L -> 7T+7T-) -2 v / 2er/+_A 



tree 



q 2 + 2pi ■ q q 2 + 2p\ ■ q 

E 2 = -^i(pi->p 2 ) • (8) 

From the dominant AI = 1/2 weak Hamiltonian, the 1 = amplitude for K$ — > 7T7t 
decays is given by 

< cc = V2G s F(Ml - Ml) , (9) 

with G s = 9.1 • 10~ 6 GeV~ 2 and F = F n = 92 A MeV at tree level. The current values for 
the CP-violating quantity r/ H are |T2 



|77+_| = (2.276 ± 0.017) ■ 10" 3 

$+_ = arg r? + _ = (43.3 ± 0.5)° . (10) 
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The magnetic tree- level amplitude starts at 0(p 4 ) where it is completely given in 
terms of low-energy constants of the nonleptonic weak Lagrangian of 0(p 4 ) |T5|, WM. In 



the notation of Ref. |15|, the leading-order magnetic amplitude is written as 



M = -^(a 2 + 2a 4 ) (11) 



with dimensionless coefficients a 2 ,a 4 of order one. As shown in Ref. , the chiral 



anomaly induces positive contributions to these coefficients. If the anomaly- induced con- 
tributions were the dominant ones, we would expect a 2 + 2a 4 to be positive. 

It is easy to check for the tree-level amplitude that the CP-violating quantity 
lm[(Ei — E 2 )M*} has a definite sign that equals the sign of a 2 + 2a 4 . From (Q) and 
(0) we would therefore conclude that a 2 + 2a 4 is positive supporting the hypothesis that 
the contributions from the chiral anomaly dominate the magnetic amplitude (|TT1). 

However, this cannot be the whole story. Assuming a constant magnetic amplitude, 
one can extract this amplitude from the branching ratio for Kl — > 7r + 7r~7. Following 
Sehgal and Wanninger [|IJ (the same procedure is used in Refs. ||, [J), one obtains in 
their notation 

with \f s \ = 3.9 • 10 7 GeV denoting the absolute value of the Ks — > 7r + 7r amplitude. 
The resulting asymmetry at tree level is 

= 7.7 % , (13) 

much smaller than the experimental result (Q). Clearly, higher-order corrections are nec- 
essary to understand the asymmetry. 

4. Chiral corrections to the tree-level amplitudes were already considered in Ref. |J. 
Those authors included one-loop corrections to the absorptive parts of both electric and 
magnetic form factors. Whereas the corrections are small for the magnetic amplitude they 
are sizable for the electric form factors, mainly due to the large final state interactions 
of two pions with / = in an S- wave, as first pointed out by Sehgal and Wanninger M. 
To some extent, Elwood et al. || compensated the neglect of dispersive contributions by 
using the tree-level value of the weak coupling constant Gs also in the loop amplitude. 
On the other hand, this is partly double counting because Gs should be reduced by about 
30 % [[ll], when effects of 0{p 4 ) are included in the amplitude for — > 7r + 7r~ . 



We propose here a different procedure for the CP-violating electric amplitude. We 
first decompose the electric form factors into a Bremsstrahlung part and a direct-emission 
piece: 

E^Ef + E™ (i = l,2). (14) 
In the first part, we use the phenomenological Ks — > 7r + 7r~ amplitude, 



Ef = 2 ^ e7] " 



q 2 + 2pi ■ q 



A e< + -^A 2 e< 
a/2 



E^ = -^(pi-Pa), (15) 
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with isospin amplitudes A = 2.71 • 1(T 7 GeV, A 2 
iment. The pion-pion phase shifts 8®, 5$ at s 



Ml 



0.12 • 10~ 7 GeV taken from exper- 
are taken from a new analysis of 
pion-pion scattering combining dispersion theory (Roy equations) with chiral perturba- 

@: 

39.1° 



tion theory |T 



5° 
°o 



51 



-8.5° . (16) 

For the CP-violating direct-emission part, we have calculated the amplitude of 0(p 4 ) 
for K® — > 7r + 7r _ 7*, including both a local amplitude derived from the nonleptonic weak 
Lagrangian of 0(p A ) [14| and a loop amplitude restricted to the dominant pion loops. 
Of course, one has to project out the Bremsstrahlung part of 0(p A ) that is already 
included in Ef in (|T^). Altogether, we have 



DE 



£r loop v 

, 2er/+_G 



8 2 



3F 



E™ 



-El 



F 
DE (P 



q- [N[ A (M p ) - N[,{M p ) - 3(iV[ 6 (M p ) + N l7 )\ 

p 2 ■ q [N r u (M p ) - N[ b {M p ) - N[,{M p ) - N 17 ) 
■+P2) ■ 



(17) 
can be 



The explicit form of the direct-emission pion-loop amplitude i?f~ loops (^ 
found in Ref. [P0 |. Because we have only included the dominant pion loops there is a 
residual, but numerically unimportant scale dependence in (|T7j). We have chosen the usual 
renormalization scale /i = M p . The numerical values for the low-energy constants N[(M p ) 
(Ni 7 is scale independent) are taken from a recent analysis 0, |21| of the branching ratio 
B(Kl — > 7r + 7r~e + e~). Anticipating the numerical discussion, the direct-emission form 
factors Ef E are negligible for Acp, especially when integrated over the whole phase 
space. 

For the magnetic amplitude, we must be less ambitious for the time being. Already at 
0(p A ), we cannot claim to be able to calculate the coefficients in Eq. (0). Moreover, 77 — 77' 
mixing is known to produce a big contribution of 0(p & ) p2| that interferes destructively 



with the amplitude induced by the chiral anomaly. In addition, there are two apparently 
non-equivalent versions of implementing vector and axial- vector exchange in nonleptonic 
weak transitions [23, 23]. 

On the other hand, there exists strong experimental evidence for higher-order effects 
parametrized in terms of a p-dominated form factor 0, |25|. Including final state interac- 
tions appropriate for P-wave irn scattering, we therefore adopt the magnetic amplitude 
measured by KTeV, 



1 + 



ai/a 2 



Ml 



Mf< + 2M k e; 

with E* the total lepton energy in the kaon rest frame and with[] [[J 



g M1 = 1.35 



+0.20 
-0.17 



ai /a 2 = (-0.720 ± 0.028) GeV 2 



;is) 



(19) 



"The quantity 02 in (IS) must not be confused with the same expression in (Ul). 
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The I = J = 1 phase shift is parametrized as 



8l(s n ) = - arcsin 



(20) 



with 



,1/2 



0.001/M^ 



(21) 



a\ = 0.038/M2+ , b\ = 0.0057/M*+ 
q* = {sjA - 

The slope parameters a\,b\ are taken from Refs. [TP| and the coefficient c\ is intro- 
duced to reproduce the correct phase shift at s — M\. 

With the sign convention of ([15]), the magnetic coupling cjmi must be positive in order 
to reproduce the measured sign of the asymmetry. Unfortunately, there is at present no 
unique way to infer from (|18|) and the experimental values (|l^) the sign of the lowest-order 
combination 02+20,4 in (p+|). A recent analysis of K — > nn^ transitions by D'Ambrosio and 



Gao p6| in the framework of the vector-field representation for (axial-)vector resonances 
finds that both signs are possible depending on the coupling strength of spin-1 exchange. 
A deeper understanding of higher-order effects in nonleptonic weak interactions will be 
necessary before one could claim that the chiral anomaly is largely responsible for the 
measured sign of the asymmetry. 



5. With electric form factors given in Eqs. fllq) and fll7j ) and with the magnetic ampli- 
tude (TT8|) we obtain for the total integrated asymmetry 

Acp = 13.7 % (22) 

if all input quantities are taken at their respective mean values. 

This value agrees with previous theoretical estimates and with the experimental re- 
sults from KTeV and NA48. The main issue we want to address here is the theoretical 
uncertainty of this prediction. Electric and magnetic amplitudes are on quite a different 
footing in this respect. Except for the P-wave phase shift SKs^) that we will lump to- 
gether with the phases occurring in Ef , we cannot ascribe a meaningful theoretical error 
to the magnetic amplitude. We will therefore scale the prediction ( p2|) to the measured 
magnetic coupling g M1 and include the experimental error of the ratio ai/a 2 explicitly. 
Accounting for the measured branching ratio in the same way, we arrive at 

B(K L — > 7r + 7r e+e ) 1.35 

where the given error is due to the error of 0,1/0-2 in ( |T9"D only. Future experimental 
improvements and correlations between the measured values for B(Kl — > 7r + 7r~e + e~), 
gMi and Oi/o 2 can easily be incorporated in this formula. 

In contrast, the uncertainty of Acp related to the electric form factors is fully under 
control theoretically. For the total asymmetry under consideration here, the CP-violating 
direct-emission amplitude (0) is completely negligible and does not affect Acp to the 
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accuracy given. Of course, this is a consequence of the dominance of the Bremsstrahlung 
form factors flT5] ) at small q 2 . 

Inclusion of the 1 = 2 amplitude A 2 is itself a small effect contributing 0.25 % to Acp- 



The relatively big errors of A 2 and <5 2 due to isospin violation [27| and electromagnetic 



corrections |28| have therefore no impact on Acp- On the other hand, those corrections 
are small for Aq and 5° that do matter for the asymmetry. Assigning an uncertainty 
±0.05 • 10~ 7 GeV to A induces AA C p = 0.25 %. Likewise, the error for |^+_| given in 
(|1C|) propagates into AAcp = 0.10 %. 

Because of the dominance of Bremsstrahlung, the phases relevant for Acp enter effec- 
tively in the combination <5q — Sl(s n ) + From the recent analysis of pion-pion phase 
shifts |19| we extract an uncertainty of 0.6° for <5q — <5j (s^). Together with the error 
of displayed in (|H]), the conservative uncertainty A (<5q — ^(s^) + $+_) = 1° gives 
rise to AA C p = 0.05 %. 

Combining the uncertainties from Aq, \tj^ | and from the phases in quadrature, we 

add the resulting theoretical error of 0.3 % to the previous one due to a\ja 2 to obtain 
our final result for the total asymmetry: 

Acp = '• 63 "+°~_ 7 + ■ (13.7 ± 1.3 ± 0.3) % . (24) 
B[Kl — > 7r + 7r e + e ) 1.35 

Given the magnetic amplitude and the branching ratio B(Kl — > 7r + 7r~e + e~), the CP- 
violating asymmetry can be calculated to a relative accuracy of about 2 %. 

6. Following the authors of Refs. || |^, [/J, we consider the asymmetry also for different 
cuts in q 2 , the invariant mass squared of the leptons. For this purpose, the chiral amplitude 
to OQo 4 ) R |2T| is used to calculate the rate T(K L — > it + -n~e + e~) in the denominator of the 



asymmetry for the cuts in question. As already mentioned, the low-energy constants N£ 
are taken from a recent analysis [[21], [H]] of the total branching ratio B(Kl — > 7r + 7r~e + e _ ). 

Suppressing errors, the results are collected in Table [I]. One immediate conclusion is 
that the asymmetry is maximal for q 2 > (10 MeV) 2 . 

Compared to the previous situation for the total asymmetry with q 2 > Am 2 , there are 
now additional uncertainties that increase with the lower limit on q 2 . In the numerator of 
the asymmetry, the direct-emission form factors Ef E in flTTD become more important in 
comparison with the Bremsstrahlung amplitude. Likewise, the rate in the denominator 
becomes more sensitive to the inaccurately known weak low-energy constants N? as one 
moves to larger q 2 . 

For large q 2 , both the asymmetry and the rate decrease rapidly. For illustration, we 
therefore choose a realistic cut of q 2 > (40 MeV) 2 . The CP-conserving amplitude of 0(p 4 ) 
depends on the combination Z(fi) = A^[ 4 (/i) — A?"{ 5 (/i) — 3(A?"[ 6 (/i) — Nn) of weak low-energy 
constants. With the mean value Z(M P ) = 0.023 extracted from the total branching ratio 
2T|, about 13 % of all events (almost 16 % of those with q 2 > (2 MeV) 2 as in the KTeV 



experiment 0) satisfy q 2 > (40 MeV) 2 . 

However, these numbers are rather sensitive to the precise value of Z(M p ). In Table 
0, the dependence of the rate on the electric amplitude is exhibited for q 2 > (40 MeV) 2 . 
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Table 1: The asymmetry Acp for different cuts in q 2 , the invariant mass squared of the 
lepton pair. 



q 2 > [MeV 2 ] 


A C p in % 


entire phase space 


13.7 


2 2 


15.2 


10 2 


15.7 


20 2 


14.0 


30 2 


12.1 


40 2 


10.2 


60 2 


7.2 


80 2 


4.8 


100 2 


3.2 


120 2 


2.0 


180 2 


0.4 



About 25% of the rate is due to Z(M P ) being different from zero. This also documents 
that Z(M p ) can in principle be extracted with much higher precision from the partial rate 
than from the total branching ratio. In the sample accumulated by KTeV [J2{|, several 
hundred events are expected to satisfy q 2 > (40 MeV) 2 . 

As already mentioned, the direct-emission amplitude of 0(p 4 ) also enters the numera- 
tor of the asymmetry. In contrast to the rate where the CP-conserving electric amplitude 
of 0(p 4 ) overtakes the CP-violating Bremsstrahlung amplitude (^j of 0(p 2 ) already at 
fairly small q 2 [21]], only the CP- violating part of direct emission matters in the numer- 



ator. For the considered cut q 2 > (40 MeV) 2 this is still only a small correction to the 
leading-order Bremsstrahlung. Numerically, the correction is smaller than the uncertainty 
due to A and it will be included in the theoretical error. 

Consequently, for reasonable cuts in q 2 the asymmetry can still be predicted rather 
precisely in terms of the magnetic amplitude and the branching ratio. For q 2 > (40 MeV) 2 , 
we obtain 

^ = B(Kl _ ^It 1 ;^ (40 MW» • if (1 °- 2 ± 11 ± °- 3) % <25) 

where the first error is again due to the error of a\/a2 in (0), the second one being the 
genuine theoretical error. 

7. The electric direct-emission amplitude for the decay — > 7i + 7i^e + e~ has been 
calculated in chiral perturbation theory. The CP-conserving part is important for the 
rate because it dominates the CP-violating Bremsstrahlung amplitude already for rather 
low invariant masses of the lepton pair. With appropriate cuts in q 2 , it will be possible to 
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Table 2: Dependence of the rate T cnt {KL — > 7r + 7r e + e ) on the electric amplitude |21| for 

q 2 > (40 MeV) 2 . E^ p ' denotes the complete electric amplitude to 0(p 4 ) that depends on 
the combination Z{n) = iV[ 4 (/i) — iV[ 5 (/i) — 3(AT[ 6 (/x) — AT 17 ) of weak low-energy constants. 
The magnetic amplitude is given in flT8|). 





r cu t/r cut (0(p 4 )) 




0.60 


(p4) with Z(M p ) = 0.023 


1. 


£° (p4) with Z(M p ) = 


0.75 



extract the relevant combination of low-energy constants in the CP-conserving amplitude 
of 0(p 4 ) with good precision. 

In contrast, the CP- violating asymmetry is quite insensitive to electric direct emission. 
This allows for an accurate calculation of the asymmetry once the magnetic amplitude 
and branching ratio have been determined experimentally. For realistic cuts in q 2 , the 
asymmetry can be predicted with a relative accuracy of 2 4- 3%. 

The contribution of the chiral anomaly to the magnetic amplitude of 0(p 4 ) leads to the 
observed sign of the asymmetry. However, higher-order terms in the magnetic amplitude, 
clearly required by experiment but not reliably calculable in chiral perturbation theory, 
make this connection less conclusive. 



We thank T. Barker, B. Cox, G. D'Ambrosio, M. Jeitler, S. Ledovskoy and L.M. Sehgal 
for helpful correspondence. 
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